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I INTRODUCTION 


NASTRAN solution techniques are shown for a numerical analysis of a 
class of coupled vector flow processes described by simultaneous parabolic 
differential equations . To define one physical problem type where equations 
of this form arise, the differential equations describing the coupled transfers 
of heat and mass in mechanical equilibrium with negligible mass average velo- 
city are presented and discussed. Also shown are the equations describing 
! seepage when both electrokinetic and hydrodynamic forces occur. Based on a 
variational statement of the general problem type, the concepts of scalar 
transfer elements and parallel element systems are introduced. It is shown 
.that adoption of these concepts allows the direct use of NASTRAN* s existing 
Laplace type elements for uncoupled flow (i.e., the heat transfer elements) for 
treating multi-component coupled transfer. Sample problems are included which 
demonstrate the application of these techniques for both steady-state and 
transient problems . 
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IRREVERSIBLE THERMODYNAMICS 


Based on the linear phenomenological laws of irreversible thermodynamics 
and the general conservation laws, generalized flow equations describing the 
simultaneous transfer of vector processes may be developed. A brief develop- 
Iment for coupled energy and mass flow is presented here and the reader is 
referred to the literature [references 1 to 8] for more details and applications. 

! Consider a motionless fluid, i.e., a fluid with zero or negligible mass 
average velocity, consisting of N chemical components K and contained in a 
mathematically well-defined geometric volume V which is separated from its 
surroundings by an equally well-defined boundary surface S . Within this system 
the coupled diffusion of matter and energy takes place. Assume that no 
chemical reactions are occurring and that the effects of external forces and 
viscous phenomena may be neglected. For this nondeforming system, therefore, 
the diffusion of momentum and its associated conservation law need not be 
directly considered. The conservation laws of mass and energy are applicable. 
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Adopting the repeated indices summation convention on i, j and designat- 
ing rectangular Cartesian coordinates as x^, the conservation of energy for 

the system under consideration can be written as 


and the conservation of mass as 
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where 

C = specific heat P K 

C = mass fraction of chemical component K(C^ - — — ) 

= i 6 * 1 spacial component of the heat flux vector with respect to the 

* mass average velocity 

jK _ ^th S p ac £ a i component of the diffusion flux vector of chemical 
i component K with respect to the mass average velocity 
t = time 

T = absolute temperature 

p = total mass density 

p = mass density of chemical component K 


The conservation equations (1) and (2) are of similar form and can be 
generalized to account for distributed external sources. Thus, a general 
conservation law for the M generalized fluxes y can be written as 
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where 

C = specific mass capacity for flux y 

Y 

jf = y spacial component of the y ^ generalized flux vector with 
1 respect to the mass average velocity 

f = scalar thermodynamic potential corresponding to flux y 

Y 

Q = external source strength corresponding to flux y 
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According to the linear laws of irreversible thermodynamics, the M 
generalized fluxes result from the action of M generalized forces X 

~E 

C e ~ 1* 2, . . . M) and are proportional to these forces. For general aniso- 
tropic coupled transfer, therefore, we may write 



where 

“ J spacial component of the e generalized thermodynamic force 
If ^ ~ phenomenological coefficient expressing the effects of the i*’* 1 

T — Hi 

component of the thermodynamic force X on the i component of 
the flux 


In addition, it can be shown that the appropriate selection of fluxes and 
i forces through examination of the local rate of entropy production leads to the 
important results: 



where 



( 6 ) 


I proportionality factor depending on the choice of generalized 

I thermodynamic fluxes and forces; tj is assumed constant for the 

| range of operating conditions considered here. 

That is, the phenomenological coefficient matrix is symmetric and the general- 
[ ized thermodynamic forces are expressible as gradients of the scalar potentials. 
When written in the above form, the quantities are known as Onsager coeffi- 
cients. 

Combining equations (3), (4) and (6) results in a general set of simulta- 
neous differential equations of the parabolic type describing the coupled 
irreversible transfer phenomena 
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(y = 1,2, .. .,M) 
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For the case of isotropic transfer equations 
respectively 


(4) and (7) become. 
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Equations (9) in their most simple form may be used to describe the 
steady-state diffusion of uncoupled heat or mass in an isotropic medium. In 
their most general form, equations (7) may be utilized to describe the transient 
diffusion of general anisotropic coupled flows such as multi- component mass 
and/or heat. 

As examples of the above developments, consider first the one— dimensional 
(i = j = 1) coupled transfer of heat and mass. For this case, 


T 1 = T U X - L 21 X 

J 1 “ L ll-ll xl— 21 

T 2 - T U X - L 21 X 

J 1 ~ L 21— 11 L 21-21 


(10a) 


2 2 X 

where J* is the heat flux and J 1 the mass flux. The coupling term L 11 

expresses the effects of the concentration gradient on the heat flux while 

establishes the effects of the temperature gradient on the mass flux. The 

21 11 
Onsager coefficient is proportional to the Soret coefficient and to the 

21 11 

Dufour coefficient. Recall that 
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The equations describing coupled electrokinetic and hydrodynamic flow in 
| porous media offer a second example. In reference 9 a finite element approach 
is used to analyze the movement of fluid relative to a solid phase and the 
j migration of ionized soil particles relative to a fluid phase, both occurring 
in an electric force field. For this case, 
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= - i 


£=>1 
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(10b) 
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where is the i th component fluid discharge velocity and J ? is the i th 


L expresses the effects 

2 


component of the current density. The coefficient 

of hydrostatic potential gradients on the discharge velocities while 

expresses the effects of voltage gradient X 2i on the current. The L^ and L* 

terms express the cross coupling between voltage gradient and fluid velocity 
and between potential gradient and electrical current, respectively. 


VARIATIONAL FORMULATION 


Equation (7) is a set of simultaneous parabolic differential equations 
describing the coupled transport processes in the continuum. While a Galerkin 
approach can be employed in a straightforward manner to evolve a finite element 
statement of the problem, the approach used here is based on a variational 
formulation similar to the principle of virtual work used in mechanics, 
i Because of this similarity, it is designated as a principle of virtual general- 
| ized work and has the advantage that physical interpretation can be associated 
| with its use. 

Consider the motionless fluid within the well-defined boundary. For each 
flow quantity y , the boundary of the volume V in which the M generalized fluxes 
are transferred may be considered as consisting of three regions S , S and S 
as follows: 12 3 

1. Over the boundary the flux is prescribed. 

2. Over the boundary S the generalized convective flux is proportional 

to the thermodynamic potential differences between the boundary values 
<y Y ) s a nd the known ambient values ( ry f ^ ^ . 
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is pre- 


3. Over the boundary S 
scribed. 

Now allow a fictitious variational change r^df^ which does not violate the 

constraints of the system; i.e., r^6f is zero over and arbitrary elsewhere. 

Also, the virtual changes are assumed to occur without the passage of time 
while the various fluxes remain constant. 


^ the thermodynamic potential ( ry f^ ) g 


For this volume consider the sum of the products of the generalized fluxes 
J Y . and the variations in their conjugate thermodynamic potentials T^Sf^, that 
is, the virtual generalized work I 
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Integration of equation (11) over the total boundary area S and application 
of Gauss' Theorem yields 
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where V-j = i component of the unit vector along the outer normal to S. 
Since the general conservation of the various transfer quantities must apply, 
equation (7) is substituted for the first integral on the right-hand side of 
equation (12), yielding 
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= 0 (13) 


Note that the surface integral appearing in equation (13) extends only over that 
portion of the boundary for which <5fy is not prescribed and that the scalar 


70 


boundary fluxes per unit surface area 
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have been employed in writing this equation. In addition, it should be noted 
that the generalized convective flux can be expressed in terms of the thermo- 
j dynamic potential differences in the form 
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where 

j hj = proportionality factor expressing the contribution from the e th 

thermodynamic potential difference to the flux J Y at the boundary 
S 2 

The characteristics of the variational statement (13) are such that its 
Euler equations and boundary conditions are those governing the case of coupled 
vector transfer considered here when the expressions (4) or (8) are employed. 

As such, it may be viewed as an alternate expression of the mathematical problem 
posed by equations (7) or (9) and the associated boundary conditions discussed 
above. These results indicate the similarities between the variational prin- 
ciple for the coupled irreversible processes and the principle of virtual work 
for mechanical systems. From this presentation the analogies between the 
variation in thermodynamic potentials (or state variables) and mechanical 
deformations, between surface flux and mechanical surface forces, and between 
thermodynamic sources/sinks and mechanical body forces are clearly indicated. 


GENERAL FINITE ELEMENT FORMULATION 


For irreversible coupled flows, the thermodynamic potentials at the grid 
points of the element are taken as the degrees of freedom. Element properties 
are determined using the principle of virtual generalized work together with 
functions expressing the state of the thermodynamic potentials within the ele- 
ment in terms of their grid point values . With this information the gradients 
of the potentials and alternately the various fluxes within the element may be 
uniquely defined in terms of the grid point values of the potentials. Assembly 
of the element equations to describe the complete system is accomplished with 
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the aid of the variational principle. This assembly is accomplished by 
requiring that the thermodynamic potentials be compatible between elements and 
that flow conservation apply at the grid points. In essence, only compatible 
states are considered and the true values of the potential are those for which 
flow conservation holds, a condition which is enforced by the principle of 
virtual generalized work. 

Details of the general finite element formulation are presented in refer- 
ence 8 and only the form of the resulting matrix equations of motion describing 
the coupled processes are presented here as follows : 


dP 

[d] {_^L } + ([ K ] + [K c J){P a } 


<y + <V + { V 


(16) 


where 

(p } = values of the thermodynamic potentials at the a grid points of 

a the assembled model 

[D] = generalized capacitance matrix for the assembled model 

[K] = generalized conduction matrix for the assembled model 

[K ] = generalized boundary convection matrix for the assembled model 

C 

{ Fq } = generalized source vector for the assembled model 

{Ft,} = generalized boundary flux vector for the assembled model 

{F } = generalized boundary convection vector for the assembled model 

Note that {P } is a column vector of length equal to the number of grid points 
ot 

of the assembled model times the number of transport quantities (i.e., a x M) . 


NASTRAN SOLUTIONS 


The previous developments are general in nature and can be applied for 
the finite element description of coupled anisotropic transfer in one, two and 
three dimensions • Once suitable interpolation functions have been chosen, the 
finite element matrices can be obtained and assembled to form the set of equa- 
tions (16) describing the idealized system using well-known procedures (refs. 8 
and 9). In what follows, the reader will be assumed familiar with these pro- 
cedures and the discussion will center on special techniques which may be 
employed to permit the direct use of NASTRAN for the analysis of M coupled 
transfer processes. 


One-Dimensional Element Formulation 

A one-dimensional element suitable for use in the analysis of two coupled 
flows is presented below and shown to be equivalent to a combined system 


?2 


of parallel one-dimensional Laplace elements, representing uncoupled trans- 
fer, connected by "scalar transfer elements," representing coupling. With 
this information the extension to the case of M coupled flows is discussed. 

Consider a one-dimensional finite element of constant cross-sectional area 
A connecting grid points £ and 0 on the x axis. Assume a linear variation in 
the two thermodynamic potential variables p..(x,t) and p (x, t) along the axis 
of the element. Thus, 2 


P-^ = = a + bx 


?2 = ^ 2^2 = c + dx 


where a(t), b(t), c(t) and d(t) are generalized coordinates which can be 
related to the four unknown grid point values of the thermodynamic potentials, 

£ Q £" Q 

i.e., p^, p^, p^, » us: >- n 8 equation (17). Evaluating the generalized coordi- 

nates in terms of the grid point potential values results in the matrix equation 
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" value of the thermodynamic potential p^ at grid point etc. 


% = length of element (£ = x - x r ) 

0 5 

x^,Xq = coordinates of element end points 


For two coupled flows, the Onsager coefficient matrix can be written as 



(19) 


Using equations (18) and (19), a generalized conduction matrix for the element 
is found in reference 8 as 
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The one-dimensional element for the analysis of two coupled processes has 
two grid points and two degrees of freedom at each of these points. This 
representation results in the four by four element matrix (20) which may be 
partitioned into the form 
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j 0 

Examination of these indicates that [K-q] represents the two by two element 

fiatrix describing the uncoupled transfer of through the thermodynamic 
potential gradient in p^ * nif^ Similarly, [K 22 ] is the two by two matrix 

lescribing the uncoupled transfer of through the gradient in P 2 = ri^f^. 

/hile 1 is the two by two matrix representing coupling between the two 


»rocesses. At this point, it is useful to introduce the concept of a 
'scalar transfer element" which may connect any two degrees of freedom and 
Ls defined by the two by two generalized conductivity matrix 




Scalar transfer element 
conductivity matrix 



(25) 


/here 


k = conductivity of the scalar transfer element 


i?ith the above observations and use of the scalar transfer element concept, it 
:an be shown that a "parallel element system" composed of two one-dimensional 
Laplace elements (i.e., NAS TRAN ROD heat transfer elements) connected by four 
scalar transfer elements (i.e., SCALAR CONDUCTION elements), see Figure 1, can 
!>e utilized to represent the single coupled transfer element of equation (20) . 

This new model has four grid points with one degree of freedom at each 
point. Grid points £ and 6 are connected by ROD 1 of area A and length 

l x q ~ x g. The degree of freedom at grid point £ is p-^ and at grid point 0 

Q 

is p-^. The element properties of ROD 1 are those of equation (22) and repre- 
sent the uncoupled transfer of J^. ROD 2 is also of length £ and area A (since 
[the geometric location of £ * is identical to £ and 0 ' identical to 0 ) and repre- 
sents the uncoupled transfer of between p| and P 2 at grid points £' and 0' as 

given by equation (23). The four scalar transfer elements represent the effects 
of coupling as described by equation (24) . The conductivities of these are 
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Generalized conductivity of _ A 
SCALAR 1 transfer element £ 12 


Generalized conductivity of A 

SCALAR 2 transfer element £ 12 


(26) 


! 

Several advantages are gained using this assembly. One advantage is that 
Iphysical insight is gained into this complex transfer problem. Another is that 
existing Laplace finite elements, like those in NASTRAN for the analysis of 
uncoupled heat transfer processes, can be utilized for the parallel models of 
the primary flows and then interconnected by the scalar transfer elements. A 
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third advantage to be gained is that a simple and practical method to automate , 
in the computer for treating an arbitrary number of M coupled processes is to ! 
model the primary uncoupled flows by M parallel systems using existing NASTRAN 
elements and to then automatically interconnect these with scalar transfer i 
elements representing the coupling. 

Two-Dimensional Element Formulation j 

The developments will now be extended to the case of coupled two-dimen- 
sional isotropic flow. The most general element for use in the investigation 
of two-dimensional flow is a triangle of thickness t . Two coupled flows will 
be considered for a triangular element lying in the xy plane and defined by the 
three grid points (£, 6, <J>) with six unknown degrees of freedom. 

Proceeding as before, assume that the potentials p^(x,y,t) and p 2 (x,y,t) 
are distributed linearly within the element by the relations ] 


p^ = a + bx + cy ; 

(28 

P 2 = d + ex + fy i 

I 

where a(t), b(t), . . . , f(t) are again generalized coordinates. Anticipating 
the results from the one- dimensional parallel element system, it is chosen here 
to express the Onsager coefficient matrix for two coupled isotropic flows as 
the sum of three matrices j 

[L e ] = [L®] + [L®] + [L®] (29j 

where 
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(32) 


rom the previous results and the form of equations (30) and (31) , it is con- 

luded that the uncoupled transfer of J and J may be included by two 

arallel Laplace triangular elements of a type available in NASTRAN. With this 

i 6 

ipproach, only the [L^] contribution of equation (32) due to coupling need be 
presented here. 

| Using equations (28) and (32) and employing standard procedures lead to the 
plement matrix [K] coupling re P resent i- n 8 between the parallel element system 

'see reference 8) , thus 
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and 


a_ = 


Ve “ x e y <j, 


y 0, ■ y e 


A 
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C £ = X 6 - X 

v (etc., in cyclic order 

5 , 0 , *) 

area of the triangular element 
coordinates of element corner points, etc. 


With the above results a parallel element system composed of two NASTRAN ' 

heat transfer triangular elements and connected by scalar transfer elements can| 

be constructed, see Figure 2. Grid points £, 9, $ ate connected by TRIANGLE 1 i 

_ 0 < 
of thickness t and properties [L ] of equation (30), Grid points £ f , 0’, (J) f 

are connected by TRIANGLE 2 also of thickness t but with properties [L^] of 

equation (31). The two triangular elements are connected by nine scalar trans- 
fer elements with properties [^^ COU pi^ n g* obtained from equation (33). The 

extension to three-dimensional problems is not difficult since one need only 
utilize existing three-dimensional heat transfer elements in the program for 
the uncoupled parallel flows and utilize the variational principle and the 
scalar transfer elements to include the effects of thermodynamic coupling. 


SAMPLE PROBLEMS 


Two sample problems are presented to demonstrate aspects of the above 
developments. Parallel and scalar transfer elements are employed and solutions 
obtained using the NASTRAN program. 

Problem 1 - One-Dimensional Steady-State Process, Three Coupled Flows 

-> 1 r\ ~y O 

The one-dimensional transfer of J , J z , and J J along the x axis of a 
continuum due to gradients in the potentials p^(x), p^(x), and p^(x) is consid- 
ered. Find the steady-state distributions in the potentials. The continuum is 
of length L - 4.0 and cross-sectional area A = 1.0, and there is internal 
generation within the continuum given by 


is the generalized conductivity for the scalar transfer element connectin; 
degree of freedom £ with degree of freedom £ 1 , etc. 
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( 34 ) 


Q x (x) 

= -(4.6 + 1.2x) 

Q 2 ( x > 

= 24.4 - 1.8x 

Q 3 (x) 

= 1.0 - 12. Ox 


The boundary conditions are specified as 


Pi (o) 

= 0.0 

P x (4) = 

32.0 


P 2 (0) 

= 0.0 

P 2 W - 

-48.0 

(35) 

P 3 (0) 

= 0.0 

P 3 (4> « 

64.0 



and the Onsager relations have been determined as 
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From equations (9), (34), and (36) the equations describing the process are 
found as 
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The solution satisfying equations (37) and the boundary conditions is given by 

P 1 ( x ) = 2. Ox 2 ; P 2 ( x ) = -3. Ox 2 ; p 3 (x) = l.Ox 3 (38) 


The finite element model is shown in Figure 3. Figure 3(a) shows the fifteen 
grid points connected by a total of twelve rod elements of length £ = 1.0, i.e.. 
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three parallel systems each consisting of four unit-length rod elements and 

representing, respectively, the three primary flows J J 2 , and J • Fig- 
ure 3(b) shows a typical parallel system connecting grid points 2, 2 T , and 2" 
(i. e. , location x = 1.0) and 3, 3 T , 3" (i.e., location x = 2.0). The proper- 

ties of the rod and scalar transfer elements for this system are found using 
our previous conclusions and are presented in Table 1. 


The vector {F^} representing the effects of distributed sources Q^(x), 
Q^Cx), Q^(x) is obtained by simple lumping procedures and applied to the grid 


points as thermal loads. The boundary conditions (35) are satisfied by con- 
straining the boundary point potentials to their appropriate values using 
single point constraints (SPC) . 


Table 2 presents a comparison between analytical results (AN.) from equa- 
tion (38) and results from the finite element solution (F.E.). 

Problem 2 - Two-Dimensional Transient Process, Two Coupled Flows 

1 -> n 

The two-dimensional coupled transport of the fluxes J and J in a square 
continuum of planform dimension 4.0 and thickness t = 2.0 due to gradients in 
the potentials p^(x,y,t), p^(x,y,t) is considered. Find the transient distri- 
bution of the two potentials when there is internal generation within the 
continuum given by 



Q-^x.y.t) 

= [sin -jp- sin -p][t + t^] 
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(39) 


Q 2 (x,y,t) 

r . TTX -MVi rj _ 0 Z 

= -[sin — sm , J [t + 3t 

4 4 



The initial 

and boundary conditions are 

specified as 




P-^O.y.t) 

= 0.0 

p 2 (°>y> t ) = 

0.0 



P, (4,y,t) 

= 0.0 

p (4,y,t) = 

0.0 
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(40) 


Pi (x , 0 , t ) 

= 0.0 

P 2 (x,0,t) = 

0.0 



P 1 (x,4, t) 

= 0.0 

P 2 (x,4,t) = 

0.0 


The Onsager 

relations for 

the isotropic 

material have 

been determined 

as 
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kd the specific mass capacities as 


pC 1 = 0.61685 


pC 2 = 0.61685 


(42) 


From equations (9), (39), (41), and (42), the equations describing the 
rocess are found as 
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'he solution satisfying (43) and the boundary and initial conditions is 


P 1 (x,y,t) 


I 


P 2 (x,y,t) 



(44) 


The finite element model is shown in Figure 4. It consists of two parallel 
five by five grid point arrays connected by a total of sixty-four right triangular 
jlements, i.e., thirty- two in each array. The two parallel systems represent, 
respectively, the two primary flows J ^ and J ^ . Figure 4(b) shows a typical 
parallel system connecting grid points 7, 8, and 13 with 7', 8’, and 13'. The 
jroperties of the triangular elements and scalar transfer elements for this system 
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are found using our previous conclusions for triangular 
For the 7,8,13 system we find, using equation (33) for 
that 


and scalar elements, 
the scalar elements. 



This information is presented in Table 3. I 

The vector {Fq} representing the effects of the distributed sources 

Q 1 (x,y,t) and Q 2 (x,y,t) was obtained by simple lumping procedures. The boundarji 

and initial conditions (40) were satisfied by constraining the boundary grid 
point potentials to their appropriate values. 

Table 4 presents a representative comparison between analytical results 
(AN.) from (44) and results from the NASTRAN analysis (F.E.). Numerical 
integration was performed using 100 time steps with integration interval 
t = 0.05 and maximum error was less than five percent. 


82 


REFERENCES 


1. DeGroot, S. R. , and Mazur, P.: Non-Equilibrium Thermodynamics . North- 

Holland Pub. Co., Amsterdam, 1962. 

2. Fitts, D. D.: Nonequilibrium Thermodynamics. McGraw-Hill Book Co., Inc., 

1962. 

3. Veinik, A. I.: Thermodynamics. Israel Program for Scientific Translations, 

Jerusalem, 1964. 

4. Taylor, S. A., Gary, J. W. : Linear Equations for the Simultaneous Flow of 

Matter and Energy in a Continuous Soil System. Soil Sci. Soc. Proceedings, 
Jan. - Feb. 1964. 

5. Lykov, A. V., and Mikhailov, Yu. A.: Theory of Heat and Mass Transfer. 

Israel Program for Scientific Translations, Jerusalem, 1965. 

6. Schechter, R. S.: The Variational Method in Engineering. McGraw-Hill 

Book Co., New York, 1967. 

7. Tsoi, P. V.: Thermodynamics of Irreversible Processes and Derivation of a 

System of Differential Equations for Molecular Transfer. J. Eng. Phys., 
Vol. 9, No. 3, 1967. 

8. Mason, J. B.: Finite Element Analysis of Coupled Irreversible Vector 

Processes. Ph.D. Thesis, University of Maryland, College Park, Maryland 
August 1971. 

9. Lewis, R. W. , and Garner, R. W.: A Finite Element Solution of Coupled 

Electrokinetic and Hydrodynamic Flow in Porous Media. Int. J. 

Numerical Methods Eng., Vol. 5, 1972. 


83 



Table 1 


Element Connection Table for Parallel 
System Between Grid Points 2, 2', 2" and 3, 3’, 3 


Element 

Connecting 

A 

£ 

L 


k - I j L 

Type 

Grid Points 



ye 


£ ye 

ROD 1 

2-3 

1.0 

1.0 

L ll“ 

1.0 


ROD 2 

2'- 3' 

1.0 

1.0 

L 22 = 

4.0 


ROD 3 

2 »._ 3 «» 

1.0 

1.0 

L 33 

2.0 


Scalar 1 

2-2’ 

1.0 

1.0 

L 12 = - 

•0.1 

0.1 

Scalar 1 

2 - 2" 

1.0 

1.0 

L 13 = 

0.2 

-0.2 

Scalar 1 

2’- 2" 

1.0 

1.0 

L 23 = 

0.3 

-0.3 

Scalar 1 

3-3' 

1.0 

1.0 

L 12 = ' 

-0.1 

0.1 

Scalar 1 

3 - 3" 

1.0 

1.0 

L 13 = 

0.2 

-0.2 

Scalar 1 

3'- 3" 

1.0 

1.0 

L 23 = 

0.3 

-0.3 

Scalar 2 

2-3' 

1.0 

1.0 

L 12 = 

-0.1 

-0.1 

Scalar 2 

2'- 3 

1.0 

1.0 

L 12 = 

-0.1 

-0.1 

Scalar 2 

2 - 3" 

1.0 

1.0 

L 13 = 

0.2 

0.2 

Scalar 2 

2"- 3 

1.0 

1.0 

L 13 = 

0.2 

0.2 

Scalar 2 

2 »_ 3 .. 

1.0 

1.0 

L 23 = 

0.3 

0.3 

Scalar 2 

2"- 3' 

1.0 

1.0 

L 23 = 

0.3 

0.3 


84 


*o|> 


Table 2 

Problem 1 Results 


Location 

P 

L (x) 

P 2 

(x) 

P, 

3 (x) 

X = 

AN. 

F.E. 

AN. 

F.E. 

AN. 

F.E. 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1.0 

2.0 

2.0“ 

-3.0 

-3.0 

1.0 

1.0 + 

2.0 

8.0 

8.0“ 

-12.0 

-12.0 

8.0 

8.0 + 

3.0 

18.0 

18.0“ 

-27.0 

-27.0 

27.0 

27. 0 + 

4.0 

32.0 

32.0 

-48.0 

-48.0 

64.0 

64.0 


Table 3 

Element Connection Table for Parallel System 



Connecting Grid 

Points 

7,8,13 

and 7 ',8’ 

,13’ 


Element 

Type 

Connecting 
Grid Points 

t 

A 

L 

ye 

k = -K , . 

coupling 

TRIANGLE 

1 7-8 -13 

2.0 

0.5 

2.0 



TRIANGLE 

2 I'-S'-O' 

2.0 

0.5 

4.0 



Scalar 

7 -7’ 

2.0 

0.5 

1.0 

“ k 7 , 7 ' = 

-1.0 

Scalar 

8 -8' 

2.0 

0.5 

1.0 

“ k 8, 8 ' 

-2.0 

Scalar 

13 -13’ 

2.0 

0.5 

1.0 

“ k 13,13' 

= -1.0 

Scalar 

7 -8' 

2.0 

0.5 

1.0 

“ k 7 , 8 ' = 

1.0 

Scalar 

8 -7* 

2.0 

0.5 

1.0 

“ k 8,7* 

1.0 

Scalar 

8 -13' 

2.0 

0.5 

1.0 

“ k 8,13' 

1.0 

Scalar 

13 -8 ? 

2.0 

0.5 

1.0 

“ k 13, 8' = 

1.0 

Scalar 

7 -13’ 

2.0 

0.5 

1.0 

“ k 7, 13 ' 

0.0 

Scalar 

13 -7’ 

2.0 

0.5 

1.0 

“ k 13 , 7 * = 

0.0 



Table 4 

Problem 2 Results 


Time 


v— N 

rt 

'w' 



P 2 (x, 

-y,t) 



x = 2.0j 

. y = 2.0 

x = 3.( 

),y = 3^0 

x = 2.( 

), y = 1.0 

x = 1 

.0, y = 1.0 

AN. 

F.E. 

AN. 

F.E. 

AN. 

F.E. 

AN. 

F.E. 

t-1.0 

0.811 

0.832 

0.573 

0.588 

-0.811 

-0.839 

-0.405 

-0.420 

t=2.0 

3.242 

3.361 

2.293 

2.376 

-3.242 

-3.381 

-1.621 

-1.690 

t=3.0 

7.295 

7.596 

5.516 

5.371 

- 7.295 

- 7.629 

-3.648 

-3.814 

t=4.0 

12.969 

13.538 

9.171 

9.573 

-12.969 

-13.584 

-6.485 

-6.792 


■ ROD 


Pl- 




<_S 


: r , x SCALAR #2 

* <V. 

1 v. i 


SCALAR ® 1 — 


*V,7 i-,. . v\Y- v 

i ^ 

*#•> 


^ si 





ROD 2 


NOTE: 

Geometric Coordinntes of Grid Points ^, Q 
Correspond to those of Respectively 


Figure 1. ONE DIMENSIONAL PARALLEL ROD ELEMENT 
SYSTEM FOR TWO COUPLED FLOWS. 


p l* 


SCALAR # 1 


p 2 Q 
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Note: Geometric Coordinates of Grid Points 
0, $ Correspond to Those of 
Respective! y 


(a) Showing Typical Scalar Connections K f*' 

TRI ANGLE 1- 




<P 


\ 


N 


:v- 


S '--' 
<3 ^r r 


/ <7:1 / 


S’ -T- ■‘% S V/.'V r V V / v . ! 

«3 .0f>' / *5? «S' 

ij %?. / £? » 

£ Tfc .,.^ivv / «' 5i 


/' 


SCALAR ELEMENTS 


M- i 


8 ,*^M*** 


^.5 


v 

f, ft 
iO'-' !.’ 


•? g 



TRIANGLE 2 



(b) Showing Al! Staler Transfer Elements. 

Figure 2. TWO DIMENSIONAL PARALLEL TRIANGULAR ELEMENT 
SYSTEM FOR TWO COUPLED FLOWS. 
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(a) Assembled Mode;! Showing Grid Points and Parallel Triangles. 
(Viscous Dampers Not Shown) 



SCALAR ELEMENTS 
NOTE: 

Geometric Location of Points 1 Through 25 
Corresponds to Those of 1* Through *25, 
Respect ively 


(b) Showing Typical Parallel Element System. 


Figure 4. TWO DIMENSIONAL, TWO COUPLED FLOWS. 
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